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Radio and gamma-ray emission from active galactic nuclei (AGN) are thought to share a common 
origin, related to the ejection phenomena in the vicinity of supermassive black holes. Thanks to its 
sensitivity, surveying capability, and broad energy range, the Large Area Telescope (LAT) onboard 
the Fermi Gamma-ray Space Telescope has permitted us to discover and characterize a huge number 
of extragalactic 7-ray sources. Similarly to what was found by EGRET, these sources are typically 
associated with blazars, characterized by significant radio emission and flat spectrum. The radio 
luminosity distribution is extended over 7 orders of magnitudes, with flat spectrum radio quasars 
clustered at higher powers and BL Lacs more scattered; the average spectral index is consistent with 
a = 0, although a few remarkable sources have a > 0.5 {S{v) oc A comparison of the radio 

flux density and the gamma-ray photon flux is presented, although claims on its signiflcance require 
a detailed discussion and Monte Carlo simulations which will be presented in a future paper. 



I. INTRODUCTION 

Around 10% of active galactic nuclei (AGN) are 
strong sources of radio emission. This includes ra- 
dio galaxies, radio quasars (flat or steep spectrum), 
and BL Lac type objects. All these sources are gener- 
ally referred to as radio loud (RL) AGN, whereas flat 
spectrum radio quasars (FSRQ) and BL Lac objects 
are collectively known as blazars. Interestingly, the 
vast majority of identified ex trag alactic sources in the 
third EGRET catalog [3EG, belong to the blazar 
class. 

RL AGN, and blazars in particular, are the most 
numerous and most luminous class of extragalactic 
gamma-ray sources. That these sources are bright in 
both gamma-ray and radio suggests a connection be- 
tween the emission processes in the two energy bands. 
Radio emission from blazars (and RL AGN in gen- 
eral) is generally accepted to be synchrotron radiation 
emitted by relativistic electrons, while the physical 
processes responsible for the 7-ray emission are much 
less well constrained. The presence of non-thermal 
synchrotron emission implies the existence of a pop- 
ulation of relativistic electrons. In the presence of 
low energy seed photons and relativistic beaming. In- 
verse Compton (IC) up-scattering of the photons is 
frequently invoked to explain the 7-ray emission. 

Indeed, in the well known blazar sequence [lol . [l^ 
[l^ . it is proposed that the synchrotron and inverse 



Compton mechanisms give rise to a connection be- 
tween the radio luminosity and the peak frequen- 
cies and relative intensities of the characteristic two- 
humped spectral energy distribution (SED) of blazars. 
On the other hand, evidence for a direct correlation 
between radio and 7-ray flux density or luminosity 
has been widely debated since the EGRET era [e.g. 

but not conclusively demonstrated, when 
all the relevant biases and selection effects are consid- 
ered [e.g. [13. 

The Large Area Telescope (LAT) onboard the 
Fermi Gamma-ray Space Telescope, with its large field 
of view and unprecedented sensitivity, is now putting 
us in the condition of a better understanding of the 
extragalactic 7-ray source population. In anticipation 
of the launch of Fermi, large projects in the radio 
band have been undertaken [e.g. [Il [H HO, HI . The 
results of these projects can now be exploited to gain 
insights into the radio properties of this population 
[see also [23 | and into the relation between radio and 
gamma-ray properties. 

In the present paper, we report in [jlT]the main re- 
sults obtained after three months of sky-survey oper- 
ation [LAT Bright AGN Sample, LBAS,E| and in ^ 
some analysis of an expanded source set based on a 
preliminary version of the 1 year catalog. Conclusions 
and future plans are given in mVl We use a ACDM 
cosmology with h = 0.71, = 0.27 and flA — 0.73, 
where the Hubble constant Hq = lOO/i km s~^ Mpc^^ 
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FIG. 1: High vs. low frequency radio flux density for LB AS 
sources; the dashed hue corresponds to Slow ~ Shigh. 



and define radio spectral indices a such that the flux 
density S{v) oc i/"". 



II. RADIO/ GAMMA- RAY CONNECTION IN 
THE LBAS 



In the first three months of sky-survey operation, 
the LAT has revealed 125 non-pulsar bright sources 
at |6| > 10° with test statistic greater than 100 
{TS > 100, corresponding to a detection signifi- 
cance of ^ lOcr); 106 of these sources have high- 
confidence associations with known AGN, 10 have 
lower-confidence associations, and only 9 remain unas- 
sociated. Therefore, considering both high- and lower- 
confidence associations, the fraction of high-latitude 
bright 7-ray sources associated with radio loud AGN 
is as high as 93%. While EGRET results suggested 
that RL AGN would dominate this population, this 
degree of dominance is nonetheless amazing. More 
than 50% (96/181) of the high-latitude > 10°) 
sources in the 3EG catalog were unassociated and it 
seems that the overwhelming majority of those were 
RL AGN! 

Of the 106 LBAS sources, 104 are blazars, consist- 
ing of 58 FSRQs [including one narrow-line Seyfert 1, 
J0948-h0022,'2l,'3|; 42 BL Lacs; four blazars with un- 
known classification; and two radio galaxies (Centau- 
rus A and NGC 1275). The basic radio properties can 
be obtained from the radio catalogs used for the asso- 
ciations, i.e. CRATES [H] and BZCat [2l|. CRATES 
provides 8.4 GHz VLA data for sources brighter than 
»S'4.8 GHz = 65 m Jy, typically with subarcsecond reso- 
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FIG. 2: High vs. low frequency radio luminosity for LBAS 
sources; the dashed line shows the best fit for a sample 
of radio galaxies Filled circles: FSRQs; open circles: 

BL Lacs; triangles: blazars of unknown type; stars: radio 
galaxies. 



lution. BZCat is a multi-frequency catalog giving low 
frequency radio data, typically from the NVSS. Not 
surprisingly, given the flux and spectral selection cri- 
teria of the samples, the LBAS sources are relatively 
bright (98/106 - 92% have ^8.4 > 100 mJy) and have 
flat spectral index ((a) — 0.02 ± 0.27), which are typ- 
ical signatures of compact self-absorbed components. 



A. Low/high radio frequency comparison 

In Fig. [T] and Fig. [2l we show the comparison be- 
tween the low and high frequency radio flux density 
and luminosity for the LBAS sources. The flux density 
plane (Fig. [T]) clearly reflects the flatness of the spec- 
tral index, with all the points (but one) close to the 
•S'low = •S'high fine; the outlying point represents Gen A, 
the only source with a significant amount of extended 
radio emission dominating the low frequency flux den- 
sity. In the luminosity plot (Fig. [5]), the common 
dependence on the square of the distance stretches 
the distribution even more. Here, the interesting re- 
sult is the relative distribution of the LBAS sources 
with respect to a complete sample of radio galaxies 
[dashed line, from[l5|: except Gen A (which lies near 
the radio galaxy region), all cores of LBAS sources 
are more luminous when compared to radio galaxies 
of the same extended luminosity, which can be ac- 
counted for by a large Doppler factor for the radio 
cores of 7-ray sources. Moreover, the low frequency 
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luminosity (generally taken as a measure of the ex- 
tended, optically thin emission) is contaminated by 
the nuclear emission for LB AS sources, so that the 
actual extended emission luminosity would be lower 
than the observed P\ow (corresponding to a shift of 
the points in the direction of the parallel arrows in 
Fig. [2|). The two main consequences are as follows: 
(1) the Doppler beaming effect can actually be larger 
than apparent from the figure and (2) 7-ray sources 
can be found even in blazars with extended emission 
as low as - 10^3 WEz-^. 



B. Radio luminosity 

In terms of the radio luminosity L,- = vL{v) (calcu- 
lated aX V = 8.4 GHz), the sources in the present sam- 
ple with a measured redshift span the range lO"^^'^ 
erg s^^ < Lr < 10^^-3 erg s~^. BL Lacs and FS- 
RQs are not uniformly distributed in this interval: 
the former span a broad range of radio luminosi- 
ties (logLr.BLLac [ergs^^] = 42.8 ± 1.2) while the 
latter are more clustered at high radio luminosity 
(logLr.FSRQ [ergs~^] = 44.4 ± 0.6). Blazars of un- 
known type have low-S/N optical spectra, so a redshift 
is generally not available and their radio luminosities 
are not determined. Of the two radio galaxies associ- 
ated with objects in the LBAS, NGC 1275 is similar 
to the BL Lacs {L^ — 10*^-^ erg s^^), while Gen A lies 
at the very low end of the radio power distribution 
(Lr = lO^^-i erg s-i). 

C. Gamma-ray vs radio plane 

By combining data in the radio archives and the 
LAT measurements, it is possible to compare the 
properties of the LBAS sources at low and high en- 
ergy. The relevant plots are shown in Fig.[3]and Fig.|4] 
and concern the peak 7-ray flux vs. the radio flux den- 
sity and the radio luminosity vs. 7-ray spectral index 
planes, respectively. These plots are also discussed in 
detail in Abdo et al. [ij. 

The calculation of a simple Spearman's rank corre- 
lation coefficient p for the distribution in the flux-flux 
plane (Fig. [3]) of the 106 objects yields p = 0.42. How- 
ever, the actual significance of this result can not be 
assessed without Monte Carlo simulations, since bi- 
ases and/or selection effects modify the distribution 
of p from the ideal case. Moreover, when FSRQs and 
BL Lacs are considered separately, quite different re- 
sults are obtained (pfsrq = 0.19, pBLLac = 0.49), so 
that it is not possible to claim significant correlations 
at this stage. Furthermore, the radio and 7-ray data 
are not simultaneous. The importance of this effect 
needs to be taken into account, and the only method 
to address its impact is to apply similar tests to a 
simultaneous sample, such as the one provided from 
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FIG. 3: Peak gamma-ray flux vs. radio flux density at 
8.4 GHz for LBAS sources; the dashed lines show the 
CRATES flux density limit and the typical LAT detection 
threshold in three months. Filled circles: FSRQs; open 
circles: BL Lacs; triangles: blazars of unknown type; stars: 
radio galaxies. 

the F-GAMMA (OVRO 15 GHz, Effelsberg & IRAM 
30m 2-230 GHz) monitoring programs. First results 
using these data indicate significant correlations [23 . 
Fuhrmann et ah, Richards et ah, in prep.]. 

The other relevant plot for the comparison of ra- 
dio and 7-ray properties is the radio luminosity vs. 
7-ray spectral index plane (see Fig. HJ. The broad 
LAT energy range permits to readily reveal the sepa- 
ration between BL Lacs and FSRQs, with FSRQs at 
largest Lr and softer indices and BL Lacs at lower Lr 
and harder indices. As far as the two radio galaxies 
are concerned, NGC 1275 is similar to BL Lacs, while 
Gen A is well displaced, having a much softer 7-ray 
spectral index than other low-power radio sources. 



III. RADIO/GAMMA- RAY CONNECTION 
WITH EXPANDED SOURCE SET 

With the continuation of the LAT operation in sky- 
survey mode, the number of detected blazars is bound 
to increase. In the l^ear catalog under development 
from the LAT team 29], more than 1000 sources have 
been detected and characterized @. In particular, 
the large field of view of the LAT has allowed us to 
reach weak gamma-ray sources all over the sky and 
the search for their possible radio counterparts has 
required a huge amount of work [see e.g. [l8|. Indeed, 
while the LBAS results were based on sources with 
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FIG. 4: Gamma-ray photon index vs. radio luminosity 
for LBAS sources. Filled circles: FSRQs; open circles: 
BL Lacs; triangles: blazars of unknown type; stars: radio 
galaxies. The object in the top left corner is the gamma- 
ray source associated to CenA. 



TS > 100 in three months of survey, the 1-year pre- 
Uminary results presented here are based on sources 
with TS > 25 (significance about 5a). As long as 
there is a significant overlap with the LBAS, the to- 
tal number of sources available for analysis becomes 
several times as large, and the explored space of pa- 
rameters is also larger. Thanks to the great sensitivity 
and broad energy range of the LAT, Fermi has already 
been successful in revealing faint 7-ray sources and in 
the characterization of their diverse photon indices. 

It is therefore interesting to have a preliminary look 
at the topics discussed in the previous section us- 
ing the 1-year database. The total radio luminosity 
range of the associated sources spans seven orders of 
magnitude; the overall distribution remains similar 
to that found for the LBAS sources. Also the dif- 
ferent distributions of FSRQs and BL Lacs are con- 
firmed, with log FSRQ [ergs~^] = 44.2 ± 0.7 and 
logLr.BLLac [ergs"!] = 42.2 ± 1.2 (Fig.EJ. In partic- 
ular, the BL Lacs remain spread over a wider interval 
of radio luminosities, even with a hint of bimodality 
(note however that the counts are smaller for the BL 
Lacs since their redshift is generally more difficult to 
determine) . 

As far as the spectral properties are concerned, we 
show in Fig. [5] the distribution of the associations in 
the high vs. low frequency radio flux density and in 
Fig. [7] the corresponding histogram of the spectral in- 
dex distribution. It is readily apparent from both 
plots that the vast majority of sources have a flat spec- 



FIG. 5: Preliminary histogram of radio luminosity distri- 
bution for AGN associated with LAT sources in the 1-yr 
catalog. 



trum ((a) = 0.06 ± 0.23), including sources as weak 
as a few tens of mJy. However, a small but significant 
tail of steeper spectrum sources is also present, which 
are described in more detail in dedicated works, such 
as M87 [l,[lH and CenA 0, as well as in the review 
by Cavazzuti et al. [11 . 

Finally, the radio flux density vs. mean 7-ray flux 
plane becomes more populated, as shown in Fig. [S] 
The Spearman's rank correlation coefficient for the 
full sample changes to p = 0.57. Of course, the 
remarks noted above still apply: different source 
classes sample different regions of the plane, data are 
not simultaneous, and observational biases could be 
present. A more rigorous analysis is therefore neces- 
sary to claim the significance and to discuss the physi- 
cal implication of the observed distribution. This will 
be the subject of a dedicated paper [Abdo et al., in 
preparation] . 



IV. CONCLUSIONS AND OUTLOOK 

Thanks to its unique capabilities, Fermi has re- 
vealed that the bright gamma-ray extragalactic sky 
is dominated by radio loud AGN, and blazars in par- 
ticular, even more than it could be probed during the 
EGRET era. Therefore, the 7-ray data from the LAT, 
as well as the radio data from large ongoing projects, 
represent a great resource to constrain the physical 
processes related to ejection phenomena in AGN. The 
preliminary analysis presented in this work certainly 
needs a more detailed follow-up, although the fraction 
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FIG. 6: Preliminary distribution of AGN associated with 
1-yr LAT sources in the high vs. low frequency radio flux 
density; the dashed lines correspond to a = ±0.5. Red cir- 
cles: FSRQs; blue: BL Lacs; black: radio galaxies; green: 
other AGN. 



and type of associated sources, the ranges of their core 
and extended radio luminosity, and the distribution of 
the sources in the various radio~7-ray planes look al- 
ready quite interesting. 

In particular, the discussion of variability and the 
use of simultaneous data, the comparison to radio data 
at higher frequencies [H and/or resolution 



251, and 



a deeper analysis with Monte Carlo simulations are 
necessary before significant claims can be made. Fi- 
nally, it is worthwhile to remind that as we discover 
weaker 7-ray sources, we will need to broaden our con- 
straints on possible counterparts, including radio-faint 
BL Lacs (HBLs), steep spectrum radio sources, and 
Seyfert galaxies. 
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the gamma-ray flux vs. radio flux density plane. Red cir- 
cles: FSRQs; blue: BL Lacs; black: radio galaxies; green: 
other AGN. 
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